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1 Introduction 



Studying the M^-boson physics is an important way of testing the Standard Model (SM) 
and searching for "new physics" . It can be done in both the electron-positron and hadron 
colliders. In e^e~ collisions, the main source of W bosons is the process of W^-pair 
production. This process was one of the most important subjects of the LEP2 experiments 
at CERN, run in the years 1996-2000, see e.g. Refs. PIE]- It also belongs to the main 
topics of a research programme of future linear colliders (LC), see e.g. Ref. |3|. In this 
process one can measure precisely the W^-boson mass and width as well as non-abelian 
triple and quartic gauge-boson couplings^. In hadron colliders (proton-proton or proton- 
antiproton), the main source of W bosons is the process of single- ly production. The 
most precise measurements of the W^-boson mass and width in hadron colliders come 
from this process, see e.g. Ref. It can also be used to extract parton distribution 
functions (PDFs) and to measure parton luminosities [5]. 

Among radiative corrections that affect the l^-boson observables considerably is the 
photon radiation in leptonic W decays. It distorts Pi^-invariant-mass distributions re- 
constructed from VT-decay products in e^e~ experiments [HIIZI or ly-transverse-mass dis- 
tributions obtained in hadron-coUider experiments [Hj. These distortions are strongly 
acceptance-dependent, see e.g. Refs. [SllZj. This radiation also affects lepton pseudorapid- 
ity distribution, which is the main tool for the PDFs and parton luminosities measure- 
ments in the hadron colliders. Therefore, precise theoretical predictions for the photon 
radiation in the leptonic W decays is of great importance for both types of high-energy 
particle colliders. In order to be fully applicable in a realistic experimental situation, such 
predictions have to be provided in terms of a Monte Carlo event generator (MCEG). 

The 0{a) electroweak (EW) radiative corrections in the on-shell W decays were cal- 
culated analytically a long time ago by several authors [01110111111121 • In the case of the 
W^-pair production in the e^e~ colliders, the calculations of the 0{a) EW corrections in a 
double-pole approximation (DPA) were done in Refs. ^3] and Pl] . The latter were imple- 
mented in the MCEG RacoonWW [13]. For single- production in hadronic collisions, 
the respective 0{a) EW corrections were calculated in Refs. [inilHKiniEl- The MCEG 
for this process, including pure QED 0{a) corrections, was provided long ago by Berends 
and Kleiss ^H]- A two- real-photon radiation cross section in W decays was calculated in 
Ref. ^ni- On the other hand, the MC package PHOTOS [20] provides a universal tool for 
the generation of photon radiation in particle decays up to 0{a'^) in the leading-log (LL) 
approximation. It was used in the MCEG YFSWW [23] for the simulation of radiative 
W decays for the VT-pair production process in e+e" collisions. 

To date, however, none of the existing MCEG for W-boson physics included multipho- 
ton radiation in leptonic W decays through exclusive QED exponentiation. Therefore, the 
influence of higher-order radiative corrections on the VF-boson observables was difficult 
to assess^. In this paper, we provide the first calculation of the multiphoton radiation in 
leptonic W decays in the framework of the Yennie-Frautschi-Suura exclusive exponen- 

^Actually, for the quartic gauge-boson coupling an additional gauge boson, 7 or Z, is required. 
^The calculation of Ref. ^H] requires two visible photons in a detector. 
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tiation j22I- This calculation is implemented in the MCEG for single- production in 
quark-ant iquark collisions called WINHAC [23] • It is a starting point for the full MC 
program for Drell-Yan-like single-iy production at the proton-antiproton (Tevatron) and 
proton-proton (LHC) colliders. The presented calculation as well as the respective MC 
algorithm can also be implemented in the MCEG for IV-pair production in the e^e~ 
collisions, such as YFSWW. 

The paper is organized as follows. In Section 2 we provide spin amplitudes for the 
Born-level process and for the process with single-photon radiation in W decays. In 
Section 3 we discuss the YFS exponentiation in leptonic W^-boson decays. Numerical 
results are presented in Section 4. Section 5 summarizes the paper and gives some outlook. 
Finally, the appendices contain supplementary formulae. 

2 Spin amplitudes 

In the calculation of matrix elements for the process of single-iy production in hadronic 
collisions, we use the spin amplitude formalism of Ref. ^]. In this approach, spinors are 
expressed in the Weyl basis, the vector-boson polarizations in the Cartesian basis, and 
the spin amplitudes are evaluated numerically for arbitrary four-momenta and masses of 
fermions and bosons. This evaluation amounts, in practice, to multiplying 2x2 c- number 
matrices by 2-dimensional c-number vectors. We give below the general spin amplitudes 
for arbitrary fermions in the initial and in the final states, and apply them later on to the 
single-iy production in qq collisions with leptonic W decays. 

2.1 Born level 

h 
fl 

Figure 1: The Born-level Feynman diagram for single-11/ production in fermion-antifermion 
collisions. 

The Born-level Feynman diagram for single-VI^ production in fermion-antifermion col- 
lisions 

^i) + MP2, a2) ^ W^iQ, A) (1) 

is depicted in Fig. [TJ where {pi,ai) denotes the four-momentum and helicity (o"j = ±1) of 
the corresponding fermion, while [Q, A) is the four-momentum and polarization of the W- 
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boson (A = 1, 2, 3). The fermions /i and /2 are members of SU {2)l doublets with opposite 
values of the weak-isospin third component and the pair /1/2 is the SU{3)c singlet. The 
spin amplitudes for this process, in the convention of Ref . pij , read 

MPia^, a,; A) = (p,) oo^^ ^2 S (p^, e*^{Q, X),Pi)Z^^^^^ , (2) 

where e is the positron electric charge, V/^/j is the element of the weak-mixing matrix 
(the CKM matrix for quarks, the MNS matrix for leptons^), Sw = sin^iy, with 6w the 
weak-mixing (Weinberg) angle; 

uj±{p) = ^/f±Jp\; (3) 

ew{Q,X) is the ly-boson polarization vector (* denotes the c-number conjugation); and 
S{. . .) is the spinorial string function, given explicitly in Appendix A. The above spin 
amplitudes are identical for any colour singlet of the initial fermion pair /1/2. 



± 




W 



fi 

Figure 2: The Born-level Feynman diagram for VT-boson decay. 

The spin amplitudes for the Born-level VT-boson decay: 

W^iQ, A) /i(gi, ri) + Mq2, ^2), (4) 
shown diagrammatically in Fig. are given by 

A^g)(A; n, r,) = --^^^ {qi) ^rM) r2 S (q,, ew{Q, A), g2);,_., , (5) 
where ri^2 denote the helicities of the final-state fermions, and C is the colour factor 

Q _ i ■ for quarks, 
1^ 1 : for leptons. 

The above spin amplitudes can be easily translated from the vector-boson Cartesian 
basis into the helicity basis, using the following transformations: 

MU>^ = ±) = ^ [TMiX = l)-tMiX = 2)] , 

V2 (7) 



MUX = 0) = M{X = 3). 



In the following we neglect masses of neutrinos and therefore do not consider mixing in the lepton 
sector. 
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Then, the Born-level matrix element for the single- production and decay is given by 
the coherent sum of the above spin amplitudes over the TV-boson polarizations multiplied 
by the Breit-Wigner function corresponding to the W propagator: 

where 

. „2\ _ / Mw^w '■ in the fixed-width scheme, 

\ Q'^Tw/Mw '■ in the running-width scheme. 

It is known that the fixed-width and running-width schemes are connected by an appro- 
priate rescaling of the line-shape parameters, here Mw and Tw |25] . 



2.2 0{a) corrections 

The cross section for Drell-Yan-like W production in hadronic collisions is dominated by 
the resonant single-VT process. Therefore, it can be described to a good accuracy with 
the help of the leading-pole approximation (LPA) [T^ ll7j. The non-LPA contributions 
are important only for specific high- VT- invariant-mass observables (e.g. in "new physics" 
searches). In this paper we concentrate on the resonant W production; the non- resonant 
contributions will be included later on. The 0{a) EW radiative corrections to the resonant 
single-lV production and decay can be divided in a gauge-invariant way into the initial- 
state corrections (ISR), initial-final interferences (non-factorizable corrections) and the 
final-state corrections (FSR), see e.g. Refs. [inillZI- The leading ISR (mass-singular) 
QED corrections can be absorbed in the parton distribution functions, in a way similar 
to the leading QCD corrections jHHSElEj- In general, the ISR corrections have a rather 
minor effect on the single- observables at hadron colliders 01201 • The non-factorizable 
corrections are negligible in resonant VT-boson production ^21IH1- On the contrary, the 
FSR corrections affect various W observables considerably [H]. This paper is devoted to 
the FSR, and the other corrections will be considered in the future. More precisely, our 
aim here is to give a theoretical description of the QED part of the FSR corrections in 
the framework of the YFS exclusive exponentiation. 

It is known that in processes involving the H^-bosons, the electroweak corrections 
cannot be split in a gauge-invariant way into the pure-QED and pure-weak ones. However, 
one can extract some parts of photonic corrections that are gauge- independent, see e.g. 
Refs. In this paper we follow the approach of Ref. 0|, where only the infrared- 

singular and fermion- mass- logarithmic terms are extracted from the virtual 0{a) EW 
corrections and combined with the real-photon contributions. They form the so-called 
QED-like corrections. The rest of the virtual photonic corrections can be combined with 
the genuine weak-boson corrections to form the so-called weak-like corrections. Another 
solution, based on the YFS separation of the infrared (IR) QED terms, was presented 
in Ref. jTH]. It differs from the previous one by subleading (non-log) terms. It can 
also be easily implemented in our calculations. In this approach, however, the weak- 
like corrections are slightly larger numerically. Of course, when the whole 0{a) EW 
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corrections are included these two approaches are equivalent. Since in this paper we deal 
with QED-like corrections only, we have chosen the solution of Ref. which is closer 
to the full 0{a) calculation. This, however, may change in the future when also the 
weak-like corrections are included. 

The major portion of the electroweak corrections can be taken into account by using 
the so-called scheme, i.e. parametrizing the cross section by the Fermi constant 
instead of the fine-structure constant a, see e.g. Refs. |2ZIIII|- In our case, this amounts 
to the replacement 



a = - aa, = -^^^^^^ (lo) 



in the hard-process parts of the matrix elements. 



2.2.1 Virtual and real soft-photon corrections 



The virtual QED-like correction to the leptonic VT-boson decay, extracted from Ref. jH], 
reads 



a 



QED 



(M,mO = - 

TT 



, , M 
2 ( In — 

nil 



In ^ In^ — - In — 



'11^ 



M 1, M 
h - In — 

nil 2 mi 

where M is the W invariant mass (i.e. = Q^), mi is the charged-lepton mass and 
m^ a dummy photon mass (an IR regulator). After combining the virtual correction with 
the real-soft-photon contribution, one obtains the virtual -|- real-soft-photon correction (cf. 
e.g. Refs. PUni): 

,2 



a 

71 



2 In 



-l ln— f + -ln 

mi J M 2 mi 6 



(12) 



where kg is the soft-photon cut-off, i.e. the maximum energy of the soft real photon up to 
which its contribution has been integrated over. When the above correction is combined 
with the appropriate real-hard-photon contribution integrated over the remaining photon 
phase space, one obtains the total QED-like correction to the H^-boson width piTH] 



s:tot 

''qed 



a 



TT 



77 
24 



TT' 

~3 



~ -1.89 X 10" 



(13) 



which does not contain mass-logarithmic terms, in accordance with the KLN-theorem |28[ 
EH] , and is small numerically. The above formulae were obtained in the small-lepton-mass 
approximation, mi -C M, which means that the terms 0{m1 / M'^) were neglected. 



2.2.2 Real hard-photon radiation 

Here we present the scattering amplitudes for single hard-photon radiation in leptonic W- 
boson decays using the spin-amplitude formalism of Ref. [21] and the notation introduced 
in the previous subsections. For the process 

iy±(Q, A) ^ /i(gi, n) + /2(g2, r^) + l{k, /t), (14) 
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fl fl fl 

Figure 3: The Feynman diagrams for W^-boson decay including single real-photon radiation 
(in the unitary gauge). 



given by the Feynman diagrams in Fig. El we obtain the 0{a) spin amphtudes 

M}j'{X;Ti,T2,k) = UJ^ri{qi)(^r2{q2)r2 

^ , Q/2?2-e* QwQ-^*^\ , , . 

[-rrv, rv, ^:Q-)s{que^{QAU.l,^_^, 

+ 7^S{q^,e;{k,K),k,ew{Q,X),q2); (15) 
2 k ■ qi ^ ' 'ri~T2 

Qf. 

' S {qi,ew{Q,\),k,e*{k,K),q2) 



2k ■q2 



Tl,-T2 



" 2k. Q ^ + 2k. Q ^ ^' 

where e^(A;, k) is the Kth polarization vector of the photon with four- momentum k (because 
the photon is massless, k = 1,2); Qj^, Qf^ and Qw are the electric charges (in units of 
the positron charge) of the fermions /i, /2 and the ly-boson, respectively; they satisfy 
the condition: Qw = Qfi — Q/2- The spinorial functions S{. . .) are given explicitly in 
Appendix A. The QED gauge invariance for these amplitudes means that 

M^D\^'y k) = 0. (16) 

We have checked numerically that after the replacement — * A; in Eq. (fT3|) . the values of 
the spin amplitudes are consistent with zero within the double-precision accuracy. 

Then, the matrix element for single-PF production and radiative W decay can be 
obtained through 

^«(ai, a2; ri, r^, /t) = \ . ^p\^i^ ^2; A) ^^^(A; n, T2, /t), (17) 

V ^^^vk + *7m^Iv ) ^ 

where the lowest-level spin amplitude Ai^p^ for the single- production is given in Eq. 
This matrix element is a coherent convolution of non-radiative spin amplitudes for W 
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production and radiative spin amplitudes for W decay. This means that it describes the 
photon radiation in the H^-decay stage only. 

As was noticed in Ref. ^Hj, the matrix element for the single-photon radiation in 
the Drell-Yan-like W production process can be, in the fixed-width scheme, split gauge- 
invariantly into the sum of matrix elements for radiative W production convoluted with 
non-radiative W decay and non-radiative W production convoluted with radiative W 
decay. This can be achieved by exploiting the partial fraction decomposition of a product 
of ly-boson propagators arising when photon is emitted from an intermediate VT-boson 
line jlH]. This simple decomposition, however, does not work in the running- width scheme. 
In this case we have 

1 1 
'11 11 

V ' ^ V ' 

<— production decay 

where Q and Q' are the PF-boson four-momenta before and after the emission of the 
photon with four-momentum k: Q' = Q — k. The two terms in the square brackets 
correspond to the radiative production and the radiative decay, respectively, but they are 
multiphed by the factor (1 -|- iTw/Myy)'^ . So in the case of the running- width scheme, 
the partial fraction decomposition of the VT-propagator works modulo this multiplicative 
factor. However, including the running VF-boson width in the case of the photon radiation 
off the ly-line leads to a violation of the QED Ward identity, see e.g. Refs. ISHUHI]- As 
was shown in Ref. ^Dj, in order to restore the respective Ward identity it is sufficient 
to include the light-fermion-loop corrections to the WW'j vertex. In the small-fermion- 
mass approximation this amounts to multiplying the respective radiative amplitude by 
the factor 

G'fls = 1 + ^^. (19) 

When we multiply our Eq. (fTHjl by Gfls, the factor outside the square brackets on the 
r.h.s. drops out, and we obtain the decomposition of the corresponding amplitude into 
the radiative production and the radiative decay - exactly as in the fixed-width scheme. 
Therefore, our matrix element of Eq. ()17|) for single-M^ production with radiative de- 
cays is valid also in the running-width scheme. Let us finally remark that although the 
compensating factor Gfls was derived for the pure light-fermion-loop contribution to the 
ly-boson width, the respective Ward identity is satisfied for any numerical value of Tw- 
So, in particular, one may use the radiatively corrected value of the PF-width. 

It should also be noted that in order to preserve a gauge-independent definition of the 
ly-boson mass and width beyond the leading order, one should use, in both the fixed- 
and running-width schemes, the pole rather than on-shell W mass and width, see e.g. 
Refs. ^32^.33] . 



1 (18) 
1 + iTw/Mw' 
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3 The YFS exponentiation in leptonic W decays 
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Figure 4: Production of a single W in quark-antiquark collisions with multiphoton radiation 
in ly-boson decay. 

As was mentioned in the Introduction, the main purpose of this work is to provide a 
theoretical prediction for the multiphoton radiation in leptonic H^-boson decays within the 
YFS exclusive exponentiation scheme. In this paper we consider the process of single-lV 
production in hadronic collisions at the parton level, i.e. 



qiiPi) + W^{Q) l{qi) + u{q,) + 7(A;i) + . . . + 7(M, (^^ = 0, 1, . . .), (20) 



depicted diagrammatically in Fig. |3] Here we do not rely on the small-lepton-mass ap- 
proximation, i.e. the formulae below are given for arbitrary final-state lepton masses. 
The 0{a) QED YFS-exponentiated total cross section for this process reads 




oo „ 
n=0 



d^qi d^q, 




(21) 



where 



Pi 




(22) 



here. 




(23) 
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is the soft-photon radiation (eikonal) factor and 



Y{Q,qi;ks) =2a m{Q,qi;m^) + B{Q,qi;m^,ks) (24) 

the YFS form factor, where B and B are the virtual- and real-photon IR YFS functions, 
given explicitly in Appendix B for arbitrary four-momenta and masses of charged particles. 
These IR functions are regularized with the dummy photon mass m^, which cancels out 
in their sum. The real-photon function B depends also on the soft-photon energy cut- 
off kg -C EcMi which means that it was integrated analytically over the photons with 
energies < kg. The photons with energies > kg are generated exclusively with 
the help of Monte Carlo techniques. The soft cut-off kg is a dummy parameter, i.e. the 
resulting cross section does not depend on it, which can be checked both analytically 
(e.g. by differentiating Eq. ()22|1 over kg) and numerically (by evaluating the cross section 
for different values of kg)- One of the advantages of exponentiation is that kg can be 
put arbitrarily low without causing any part of the cross section to become negative - 
in contrast to fixed-order calculations. In Eq. (j^ . /^q^"* and f3^^ are the YFS non-IR 
functions, calculated perturbatively through 0{a). We present them below in the centre- 
of-mass (CM) frame of the incoming quarks, i.e. the rest frame of W, with the +z axis 
pointing in the quark qi direction. 
The function (3q^^ is given by 

f^o\pi,P2,qhqu) = /3o\pi,P2,qhqu) [i + 6'-^\Q,qi,q^)] , (25) 

where /9q°^ is related to the Born-level cross section through 

^^(0)= ^ ^ = I (26) 

2^° ^X{l,mi/M,m,/M) d^i 16s (27r)2 12 ^ I I ' ^ ^ 

with s = (gi + ^2)^ = and \{x, y, z) = + y"^ + z"^ — 2xy — 2xz — 2yz. The factor 
J2 ~ I ' I corresponds to averaging over the initial-state quark spins and colours (the 
colour contents has been extracted explicitly), and the sum ^ runs over all the initial- 
and final-state spin indices. In Eq. ()25|). the correction 

S'^'^\Q, qi, qu) = S^wiQ^ qu qv] m^) - 2a^B{Q, qr, m^) (27) 

is the 1st order non-IR correction to the /?o function, where 5'^^^ is the 0{a) EW virtual 
correction. Since in this paper we limit ourselves to the QED-like corrections, from 
Eqs. flTT|) and we have 



5SD(Q,^0 = ^(ln^ + ^). (28) 



Although Eqs. (|TT|l and were obtained in the small-lepton-mass approximation, nii <^ 
M, we have checked that the above formula remains true for arbitrary lepton mass mi < M 
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(of course, under the assumption that 5qed contains only the IR- and mass- singular 
terms) . 

The function j3^^ is the YFS non-IR function corresponding to the single-real-hard 
photon radiation. It is related to differential cross sections through 



where 



with 



dWdk^dVtk 32s (27r)5 12 



l^El-^^'f' (30) 



|2 



^/W) = — ^ (31) 

|gl|(M-A;0) + gO|A;|cos% 

the phase-space factor (coming from the phase-space integration eliminating the energy- 
momentum conservation ^'^^-'-function for single-photon radiation), where Oik = ^(Qhk); 
in the soft-photon limit A{k 0) A. The sum ^ in Eq. (|3U|) again runs over the 
initial- and final-state spin indices, this time inluding also those of the radiative photon. 
Thus, we finally have 

A'\P1,P2, qu Qu, k) = ] ^ Yl l-^^'f - ^(^' ^)/^r(Pi, P2, qi, Qu). (32) 

lbs [l-Ky 12 ^-^ ' ' 

There are several advantages in using the matrix elements of Section |21 Firstly, the 
respective spin amplitudes are derived without the assumption of the energy-momentum 
conservation. Therefore, they can be used directly in evaluations of the above YFS 13- 
functions over the multiphoton phase space, without the need to resort to any "reduction 
procedure" , which reduces the multiphoton phase space to the 0-photon phase space for [3q 
and the 1-photon phase space for see e.g. j221ini|- Secondly, since the spin amplitudes 
are obtained for massive fermions, there is no need to use any phase-space slicing or 
subtraction methods in order to separate mass singularities ^7]. Using spin amplitudes 
instead of explicit analytical formulae for the squared matrix elements may also be useful 
for some dedicated studies, such as investigation of various VF-polarization contributions, 
"new physics" searches (spin amplitudes can be easily modified to include some "new 
physics" components), etc. And, which is important in practice, the numerical evaluation 
of the matrix elements based on the above spin amplitudes is fast in terms of CPU time. 

In computing the matrix element ^ [A^^^-*! we observed a loss of numerical precision 
~ (9(0.1%) when the angle between the radiative photon and the electron (positron), 6^-^, 
was ~ (9(10~^). It turned out that most of this precision loss was coming from huge 
numerical cancellations between the terms in the universal eikonal factor of Eq. (jl5j) (the 
factor in front of the first ^-function) . We improved this by correcting the above matrix 
element according to 

5^|A^«f ^5^|A<«f + 5,,,, (33) 
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where 



(34) 

Algebraically, the two terms in the square brackets are identical. Numerically, however, 
they can differ for ultra-coUinear photon radiation, owing to huge cancellations in the 
second term leading to a loss of numerical precision. Therefore, this correction effectively 
replaces the numerically unstable part of the matrix element ^ |7V1*^^^| corresponding 
to the second term in Eq. ()34|) with the numerically safe one corresponding to the first 
term, obtained directly from the particles four-momenta. We have checked that the 
above modification is sufficient for the numerical precision of O{10~^) for Oe-y ^ 10~^ and 
of O{10~^) for the total cross section. By looking at Eq. ()32|1 one can notice that the 
part of the matrix element ^ [A^*-^-*! that is proportional to the soft-photon factor S 

exactly cancels in the calculation of the (3[^^ function. We could, therefore, perform this 
cancellation algebraically and thus avoid the above numerical problems. We, however, 
keep this term in \-M^^^ \ because apart from the YFS exponentiation we want to have 
in our program also the non-exponentiated, fixed-order 0{a) calculation. Since it is now 
calculated in the same way as the second term in Eq. it exactly cancels numerically 
in the evaluation of the f3^^ function. 

This completes our description of the cross section for process (PUI) with the 0{a) 
QED YFS exponentiation. In order to compute this cross section and generate events, 
we have developed an appropriate MC algorithm, which will be described in detail else- 
where ^5 . We will complete this paper by presenting some results of numerical tests of 
the corresponding MC program, called WINHAC. 

4 Numerical results 

We performed several numerical tests of the MC event generator WINHAC, which im- 
plements the calculations presented above. Here we discuss some of the results. We 
considered the following process: 

d + u — >W- — + (35) 

where I = e, fi, r. We have checked that the results remain unchanged when we switch to 
the corresponding process of production and decay. Our MC calculations were done 
using the scheme and the fixed-width scheme. All the results below, unless stated 
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otherwise, have been obtained for the following input parameters: 

mrf = 3 X 10~^ GeV, m„ = 6 x 10"^, Kd = 1, = 0, 

me = 0.511 X 10"^ GeV, = 0.10565836 GeV, = 1.77703 GeV, 

Mw = 80.423 GeV, Mz = 91.1882 GeV 
2 _3G^f 2as\ (36) 

J, Ur^ = l.lUUO» ^ lU ^ ' 

EcM = y/s = Mw 
4.1 General tests 



a = 137.03599976, = 1.16639 x 10"^ GeV'^ a, = 0.1185, 



We have performed several numerical tests of the MC event generator WINHAC aimed at 
checking the correctness of the implemented matrix elements as well as the corresponding 
MC algorithm. 

In order to cross-check the matrix elements presented here, we implemented in our MC 
program the matrix elements of Ref. which in the following we shall call B&K. These 
latter matrix elements were obtained in the small-lepton-mass approximation mi -C Mw', 
their precision therefore is of 0{mf/M^), which for electrons gives (9(10~^°). Since 
our spin amplitudes are obtained for massive fermions, we performed the comparisons of 
these matrix elements for electronic VT-boson decays. We did this by taking the difference 
between the corresponding MC weights on an event-by-event basis and calculating the 
average of this difference over the whole MC sample. For both the Born-level and 0{a) 
matrix elements, we reached an agreement at the level of ~ 10~®. 

Then, we performed several tests to check the MC algorithm of the program WINHAC. 
An important test of the algorithm for MC integration and event generation according to 
Eq. ()21|1 is to reproduce fixed-order calculations. The strict Born-level cross section can 
be obtained from Eq. ()21|1 by truncating the perturbation series in a at the lowest-order 
term, which amounts to 

J Hi 1u 

Within the multiphoton MC algorithm, this means calculating an appropriate weight if 
the photon number n = and setting it to zero if n > 0. The Born-level total cross 
section can be easily calculated analytically. In the small-fermion-mass approximation 
and in the fixed-width scheme it reads 



36stv {s-M^y + Miri 



tot _ 'JM I /qo^ 



In Table^we compare the results for the total Born cross section for e, /i and r in the final 
state, calculated with the MC program WINHAC with those obtained from the analytical 
formula of Eq. (jHHj) . We see a very good agreement between these two calculations for e 
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Calculation 


(T*°* [nb] 


e 




T 


Analytical 


8.8872 


8.8872 


8.8872 


WINHAC 


8.8869 (2) 


8.8873 (2) 


8.8808 (2) 



Table 1: The results for the total Born-level cross section from the MC program WINHAC compared 
with the analytical calculation in the small-fermion-mass approximation. The numbers in parentheses 
are statistical errors for the last digits. 



and /i. For r they differ by ~ 0.1%, which can be explained by the r-mass effects (they 
are not negligible as in the case of e and //). 

In a similar way, the first-order cross section can be obtained from Eq. (j2T| by trun- 
cating the perturbative series at 0{a) beyond the Born level, i.e. 



cr 



tot 



+ 



d^qi d?qy d?k 



qi - qv) Po 



(0) 



qi 



ql 



S^'\pi+P2-qi-qu-k) ^{'^ + 5^°^ 9{k'-k 



(39) 



where the first term on the r.h.s. corresponds to the Born plus virtual and real-soft-photon 
contribution, and the second term to the real-hard-photon contribution. In practice, this 
means that the first term is evaluated within the multiphoton algorithm only for n = 0, 
the second only for n = 1, otherwise the appropriate MC weights are set to zero. In 



Calculation 


5i = af7a*°* - 1 


e 




T 


WINHAC 


-1.5(3) X 10-* 


-2.2 (3) X 10-* 


-0.3 (2) X 10-* 



Table 2: The results for the 0{a) QED-like correction to the total cross section from the MC 
program WINHAC. The numbers in parentheses are statistical errors for the last digits. 

Table |21 we show the results from the program WINHAC for the pure 0{a) QED-like 
correction to the total cross section. As can be seen, the results for e and /i are in good 
agreement with the numerical value of total QED-like correction to the IT^-boson width 
as given in Eq. ()13|1 . For r we observe the difference of ~ 1.5 x 10"*, which again can be 
explained by the r-mass effects. 

In Table El we compare the results for the 0{a) hard-photon correction as a function 
of the lower photon-energy cut-off k^, i.e. 

(fco) = 4* / dE,^ X 100% , (40) 

(^1 Jko -^7 

for the centre-of-mass energy Ecm = 90 GeV, obtained from the program WINHAC and 
from the B&K MC program \TE\. The results of these two programs agree very well within 
the statistical errors. 
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e 


fi 


WINHAC 


B&K 


WINHAC 


B&K 


0.01 


19.69 (3) 


19.7 


10.11 (2) 


10.1 


0.05 


11.61(2) 


11.6 


5.92(1) 


5.9 


0.10 


8.31(2) 


8.3 


4.22(1) 


4.2 


0.15 


6.47(2) 


6.5 


3.27(1) 


3.3 


0.20 


5.23 (1) 


5.2 


2.63 (1) 


2.6 


n 'in 
u.ou 


3.61(1) 


3.6 


1.80(1) 


1.8 


0.40 


2.57(1) 


2.6 


1.27(1) 


1.3 


0.50 


1.84(1) 


1.8 


0.91(1) 


0.9 


0.60 


1.29(1) 


1.3 


0.63(1) 


0.6 


0.70 


0.86(1) 


0.9 


0.42(1) 


0.4 


0.80 


0.52(1) 


0.5 


0.25(1) 


0.2 


0.90 


0.24(1) 


0.2 


0.11 (1) 


0.1 



Table 3: The fraction of events (in %) with a photon energy greater than ko at 0{a) from the 
MC program WINHAC and from the MC program of Berends & Kleiss [TBI (denoted as B&K) for 
EcM = 90GeV. The numbers in parentheses are statistical errors for the last digits. 

As the above fixed-order results from WINHAC have been obtained in the framework 
of the YFS-type multiphoton algorithm, they make us strongly confident in the correctness 
of the corresponding MC algorithm. 



Calculation 


[nb] 


e 


/i 


r 


Fixed 0(a)-level 


8.88564(14) 


8.88539(12) 


8.88047(10) 


YFS exponentiation 


8.88390 (6) 


8.88443 (9) 


8.87859 (9) 


Sexp = (c"yFS ~ 0"*°*)/(Jq°* 


-2.0(1) X 10"^ 


-1.1 (1) X 10^^ 


-2.1 (0) X 10-^ 



Table 4: The results for the fixed-C'(a) and the YFS-exponentiated total cross section from the 
MC program WINHAC. The numbers in parentheses are statistical errors for the last digits. 

In TablelUwe give the results for the total cross section at the fixed (9(a)-level and in- 
cluding the YFS exponentiation as given in Eq. ()21|1 . The YFS-exponentiation corrections 
beyond 0{a) are ~ 10~^, i.e. of the expected size of higher-order corrections. 

4.2 Distributions 

Here we presents the results from WINHAC for some distributions at the Born level, 
with 0{a) QED-like corrections and including the YFS exponentiation. These results are 
given for two kinds of event selection: BARE - where the corresponding observables are 
obtained from bare-lepton four-momenta and no cuts are applied, and CALO - where the 
photon four-momenta are combined with the charged-lepton four-momenta if the opening 
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BARE 



CALO 




10.0 20.0 30.0 40.0 10.0 20.0 30.0 40.0 

Figure 5: Distributions of the charged-lepton energy for BARE and CALO acceptances. 

angle between their directions Z{qi,k) < 5°; such photons are discarded. No extra cuts 
are apphed. The BARE acceptance is closer to experimental event selections for muons, 
while the CALO is closer to the ones for electrons. We, however, use them for both types 
of final states. 

In Fig. we present the distributions of the electron and muon energy for the fixed 
0{a) corrections and for the YFS exponentiation. The upper plots show the absolute 
distributions for the BARE and CALO acceptances, respectively, while the lower plots 
show the relative differences between these two calculations, also for the BARE and CALO 
acceptances. At the Born level, the charged-lepton energy is fixed at Ei ^ |v^; therefore, 
the energy tails for Ei < ^a/s in the above plots are the result of the real-photon radiation 
in H^-boson decay. The YFS-exponentiation corrections beyond the fixed 0{a)-leyel are 
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BARE 



CALO 




-0.05 



-0.10 



-0.15 




0.00 



-0.02 



-0.04 



E^GeV] 
30.0 40.0 




E^[GeV] 

.1 I . . 

30.0 40.0 



10.0 20.0 30.0 40.0 10.0 20.0 

Figure 6: Distributions of the hard est- photon energy for BARE and CALO acceptances 



large for the BARE acceptance: up to ~ 15% for electrons and up to ~ 8% for muons 
- they differ for these two decay channels. For the CALO acceptance these corrections 
are smaller, up to ~ 4%, and they are almost identical for electrons and muons. This is 
because in this case the large corrections due to lepton-mass-log terms have been excluded 
by the photon-lepton recombination. As can be seen, the largest (negative) corrections 
are in the first radiative bin (i.e. the second highest one in the upper plots) and they 
change sign for low lepton energies. 

In Fig. iniwe show the distributions of the hardest photon energy for the electron and 
muon H^-decay channels. The notation is similar to that in Fig. El Again, for the BARE 
acceptance the YFS-exponentiation corrections beyond the fixed 0{a)-leve\ are large and 
different for these two channels: up to ~ 15% for electrons and up to ~ 7% for muons. 
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For CALO they are smaller, ~ 4%, and similar in the two channels. The corrections are 
largest for soft photons and decrease with the photon energy. 

Finally, in Fig. [7|we show the distributions of the cosine of the hardest photon polar 
angle with respect to the incoming quark direction; the notation is as in Fig. IHl For the 
BARE acceptance the YFS-exponentiation corrections beyond the fixed (9(a)-level are 
large and almost constant: ~ 28% for electrons and ~ 14% for muons. For CALO they 
are smaller, 6-7%, and similar in the electron and muon channels. 

As can be seen from Figs.|SHZl the YFS exponentiation affects sizeably radiative events. 
All the above distributions have been obtained for the parton-level VT-boson production 
at fixed CMS energy. In the actual proton-(anti)proton collisions the parton-parton CMS 
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energy can change, which leads to an enhancement of the FSR corrections, in particular 
those due to the YFS exponentiation. This will be investigated elsewhere jHH] . 

5 Summary and outlook 

In this paper we have presented the calculations of the YFS QED exponentiation in 
leptonic IV-boson decays. We have provided the fully massive spin amplitudes for the 
single ly-boson production and decay, including the single-real-photon radiation in W 
decays. We have obtained the numerically stable representations of the YFS form factor 
for the charged-particle decay. All this has been applied to the process of Drell-Yan-like 
W^-boson production in hadronic collisions and implemented, at the parton level, in the 
Monte Carlo event generator WINHAC 1.0. For this purpose, an efficient multiphoton MC 
algorithm has been developed. The above spin amplitudes have been cross-checked with 
the independent analytical representations of the appropriate matrix elements JH] and 
they have been found to be in very good numerical agreement. We have also performed 
several numerical tests of the implemented MC algorithm. The results of these tests make 
us confident in the correctness of this MC algorithm. 

Numerically, the YFS-exponentiation corrections beyond the fixed 0{a) calculations 
are at the level of ~ 10~^ for the total cross section, which is the result of the KLN- 
theorem. However, for some distributions they can amount to between a few and over 
20 per cent. These corrections can be significantly reduced when a calorimetric-like re- 
combination of radiative photons and charged leptons is applied. Such a treatment is 
experimentally natural for the electrons in the final state, but less obvious for the muons. 

Here we presented the calculations for the QED-like corrections in the leptonic W- 
boson decays and for the parton-level VT-production process only. We are planning to 
extend this, in the future, to the full proton-(anti)proton collisions and to include other 
0{a) electroweak corrections. The next step would be the inclusion of the NLO QCD 
effects as well as soft-gluon resummation corrections. We are also going to perform further 
tests of the program WINHAC at the 0{a) and beyond, particularly comparisons with 
independent calculations for various observables. Last, but not least, the full documen- 
tation of the MC program WINHAC [23 is in preparation (to be submitted to Computer 
Physics Communications) . There, the details of the corresponding MC algorithm will be 
given. 

In this paper, we applied the QED YFS exponentiation in leptonic 1^-boson decays 
to the single-W^ production process at hadron colliders. However, it can also be used to 
describe the photon radiation in W decays in the processes of VT-pair production at both 
hadron and electron-positron colliders. In particular, it can be rather easily implemented 
in our MC event generator YFSWW ^T] for W^W~ production in e~^e~ collisions, which 
will be necessary for the future linear colliders [3]. 
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A Spinorial string functions 



Here we provide exphcit formulae for the spinorial string functions introduced in Section|2l 
The general such function in the two- component Weyl-spinor basis reads j21I 



S{pi, ai, . . . , a„,pj)a,<7,- = xUiPi) • • • ' Xa,{Pj), 



where 

x+{p) 



1 



\p\ +p 

+ ip'^ 



x-ip) 



-p + ip 
. \P\+P^ 



(41) 



(42) 



2\p\{\p\ +p3 

are the two-component Pauli spinors corresponding to an external fermion with four- 
momentum p = {p^,p) = {p^ ,p^ ,p'^ jP"^)] for p^ = —\p\ we choose 



x+{p) 



x-{p) 



-1 




The internal part of the above string function 

[ai, . . . ,a„]" = (fii)a(^2)- 

is the product of 2x2 c-number matrices, where 

aP =F 
=F(a^ - ia^) 



in)(-l)"+la 



=p(a — la 



(43) 



(44) 



(45) 



with a = {a^, a^, a?, a^) the four-vector in the Minkowski space. 

As can be seen, the spinorial function S can be easily evaluated numerically for arbi- 
trary n. One can just compute a product of internal 2x2 matrices ((zij)a, and then multiply 
the resulting matrix by the external 2-dimensional c-number vectors x- However, the nu- 
merical evaluation is more efficient if, instead of matrix-by-matrix multiplication, one 
performs matrix-by- vector multiplication. In our computation of the function S, we start 
from multiplying the left-hand-side vector by the matrix {^i)ai and continue by mul- 
tiplying the resulting vectors by the consecutive matrices {(/ii)a until we reach the last 
matrix, {(/Ln)a- The computation is completed by performing the scalar product of the 
final vector of the above multiplication with the right-hand-side vector x- 

Three polarization vectors of a massive vector-boson with four-momentum k = {k^, k) = 
(/c°, fc^, /c^, k'^) and the mass m are, in the Cartesian basis, given by 



e'^(A;,A 


= 1) 


e'^(A;,A 


= 2) 


e'^(fc,A 


= 3) 



\k\kT 

— * 

m\k\ 



(46) 



k^ P 
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where kx = ■>/ (A;^)^ + (A;^)^ is the transverse momentum. For massless vector bosons, 
such as photons, e^(A = 3) = 0, i.e. there are only two non-zero polarizations e'^(A = 1) 
and e'^(A = 2). Helicity eigenstates can be obtained from the above polarization vectors 
through 

v2 (47j 
ef,,i{k,X = 0) = e{k,X = 3). 

B The YFS IR functions 



In Ref. we provided the general formulae for the YFS IR functions ?R.B and B for a 
pair of charged particles of arbitrary masses and four-momenta. This representation works 
very well for particle production or scattering processes; however, it becomes numerically 
unstable for charged-particle decays. Therefore, for the process 

W^{Q)^l^iq) + ii{q'), (48) 

we have to obtain different representations for these functions. A specific feature of the 
above process is that a QED-radiation dipole is stretched between a decaying particle 
{W) and a decay product (Z). For such a dipole, a four-momentum transfer between its 
constituents is non-negative: 

t={Q-q)'> 0, (49) 

in contrast to scattering processes. In such a case, the corresponding IR integrals have 
to be calculated in a slightly different way than was done in Ref. [36j for t < 0. Special 
care is also needed for the limiting case t = 0, which occurs for the two-body VT-boson 
leptonic decay when the neutrino mass is neglected. For the sake of numerical stability, 
it has to be treated separately. 

The YFS virtual- and real-photon IR functions for a pair of charged particles with the 
four- momenta {Q, q) are defined as follows ^22] 

f d'k ( 2q-k 2Q-k V 

^W, q, m^) j p _ ^2 + _ 2kq + ie P - 2kQ + le ) ' ^ ' 

where is a dummy photon mass used to regularize the IR-divergent integrals (m^ -C 
ks), while kg is the soft-photon cut-off, up to which the integration over the real-photon 
four- momenta is carried over analytically {ks ^ Q^)- The explicit analytical formulae for 
these functions are presented below. 

B.l The virtual-photon IR function 

The virtual-photon IR function reads as follows: 
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I. t = {Q-qy>0: 

,2 



2a^B{Q, g; m^) = ^ I [ iyA{Q, q)-l]ln^ + ^A, (Q, q) - i^A^ (Q, q) | , 



with 



t m t 



+ 'LA+i:_21n" + '':'"")ln' + '^ 



4 \^ m? m? J m? 

+ ^ In 77 lii(l + ,7) - 1 In C ln(l + C) + 5RLi2(-77) - UU^i-^ 



where 



u^Qq, A = ^/{v - Mm){v + Mm), = M^, = M >m, 
t = m"^ + rn^ - 2u, Mm <^ <^ + , 

_ m'^t A _ + 

2A(2A + z/ - m2) ' 

and 

Li2(y) = -/ — ln(l-a;), | arg(l - y)| < tt, 
Jo 2: 

is the Spence dilogarithm function. 
II. t^{Q-q)^ = 0: 

2aKS(g,g;m^) = -<^ U— -1 In + - 



TT \ \M^-m^ m J \ Mm 2 
In the hmit m <^ M we get 

2aUB(Q,q;m^) = -<^21n l ln-^ + ln^ -In 

m«M TT \ V"^ J M m 2 m 2 

B.2 The real-photon IR function 

For the real-photon IR function we obtain 
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I. t = {Q-qf>0: 



a 



B{Q,q;m^,ks) = - 

TT 



with 



4P 

[uA{Q,q) -1] In ^ 



—Ai{Q,Q) - —Ai{q,q) 



uA^iQ,q) }, 

















t 



[X{z^, y,, y,, y2, ya) + R{zi) , 



(60) 



(61) 
(62) 



where 



and 



+ ^^(^; yi^ y2, Vs, ydX{z; ya, Va, yi, yd , 



Yi,{z)=2Z,,{z) + -\n' 



z-Vi 



yj 



Zij{z) 



yj - yi 



X(z; a, 6, c, d) — In 



z-y, 

[z — a){z — b) 



{z — c){z — d) 



Zq 



M 

T ' 



\Q\ 

T 



yi 2T 



1 

_ 1 

ys - 



t 
V 



^1 = ^-1; 

kV 

Z/2 = 2/1 - ^ : 



y4 = 2/3 + 



(63) 



(64) 



n= ^(^cu^ -t){S^ -t), S = M -m, uj = M + m, 

while z/ and ^(Q, q) are as given in the previous subsection. We have checked 
that this analytical representation is numerically stable for t > 10^^'^GeV^, when 
computed in any Lorentz frame, which is neither W nor / rest frame. In these frames 
we need an explicit analytical formula for A4(p,p) in the limit /3 — > 0. It reads 



p={mfl) 



m 



2 ■ 



(65) 
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In the W rest frame, i.e. for Q = (M, 0), the function A^^Q, q) can be simphfied to 

get 



MQ,q)-- 



1 



2Mq 



M-E + qE + q 2q{M - E + q) ^ E + q 
In — — in — — 2 in — in 



M -E-q M 



+23fJLi2 



M -E-q 



23?Li, 



Mm 

E + q 
M 

M -E + q 



M 



(66) 



-KLis 



. fM{E + q)- 



in 



2Mq 



3?Li2 



2q 

M{E - q) - 
-2Mq 



where E = q^ , q = \q\. 
1. t^{Q-qf = 0: 

For t = the functions A{Q^ q), A4^{Q, Q). Ai(g, q) can remain the same as for t > 0, 
but we need a new, numerically stable, representation for the function A4{Q,q). It 
can be cast in the form 



In — — In 
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+ U{zi) - U{zo) 



(67) 



where 



U{z) = lln' 
+ 23fiLi2 



{z - yi){z - y2) 



y2 - yi 
z-yi 



+ 23?Li2 



+ In I2; - yi\ In 

ys - 1/2 



\z - yi\ 
{z - 2/2)^ 



z-y2 



(68) 



yi = -r, y2 = yi- tt-t? , 1/3 = -yi + 



2m^ 



A ' 2A2 ' ' 



In the W rest frame we get 



MQ,q) 



--{M,6) 



M2 



W h L12 

m 



(69) 



Then, in the small- lepton-mass limit, m -C M, we obtain a simple expression for 
the function B in the W rest frame: 



2aB{Q,q;mj,ks) 



a 



M 



m 



2k 



M , M 



TT 



- <^ 2 In 1 In — - In^ — + In — + 1 \. 



m m 



(70) 
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After combining this with the virtual-photon function of Eq. (jSH)), we obtain a simple 
expression for the YFS form factor in the W rest frame: 

YiQ,q;Q = 2 In— -1 ln-^ + -ln— ----L 71 

As can be seen see explicitly, it is free of the IR singularity as well as of the Sudakov 
double-logarithms. 

References 

[1] G. Altarelli, T. Sjostrand, and F. Zwirner, editors, Physics at LEP2 (CERN 96-01, 
Geneva, 1996), 2 vols. 

[2] S. Jadach, G. Passarino, and R. Pittau, editors. Reports of the Working Groups on 
Precision Calculations for LEP2 Physics (CERN 2000-009, Geneva, 2000). 

[3] ECFA/DESY LC Physics Working Group, J. A. Aguilar-Saavedra et al, TESLA 
Technical Design Report Part III: Physics at an e^e~ Linear Collider, SLAC- 
REPRINT-200 1-002, DESY-01-011, DESY-2001-011, DESY-Ol-OllC, DESY-2001- 
01 IC DESY-TESLA-2 001-23, DESY-TESLA-FEL-2001-05, ECFA-2001-209, Mar. 
2001, 1h^p-ph /0106315l 

[4] G. Altarelli and M. Mangano, editors. Proceedings of the Workshop on Standard 
Model Physics (and More) at the LHC (CERN 2000-004, Geneva, 2000). 

[5] M. D ittmar, F. Pauss, and D. Zurcher, Phys. Rev. D56, 7284 (1997), 
|h ^x/9 705004 

[6] W. Beenakker, F. A. Berends, and A. P. Chapovsky, Phys. Lett. B435, 233 (1998), 
hep-ph/9805327 , 

[7] S. Jadach, W. Placzek, M. Skrzypek, B. F. L. Ward, and Z. Was, Phys. Rev. D61, 
113010 (2000). 

[8] U. Baur, S. Keller, and D. Wackeroth, Phys. Rev. D59, 013002 (1998). 

[9] W. J. Marciano and A. Sirlin, Phys. Rev. D8, 3612 (1973). 
[10] J. Fleischer and F. Jegerlehner, Z. Phys. C26, 629 (1985). 
[11] D. Bardin, S. Riemann, and T. Riemann, Z. Phys. C32, 121 (1986). 
[12] A. Denner and T. Sack, Z. Phys. C46, 653 (1990). 

[13] W. Beenakker, F. A. Berends, and A. P. Chapovsky, Nucl. Phys. B548, 3 (1999), 
,hep-ph/9811481, 

24 



[14] A. Denner, S. Dittmaier, M. Roth, and D. Wackeroth, Nucl. Phys. B587, 67 (2000), 
|hep-ph/0006307 



[15: 

[16 

[17 

[18' 
[19 
[20 

[21 

[22 
[23 

[24 
[25 
[26 
[27 
[28 
[29 
[30 
[31 
[32 
[33 



D. Wackeroth and W. Holhk, Phys. Rev. D55, 6788 (1997). 

V. A. Zykunov, Lowest order electroweak radiative corrections to the single W pro- 
duction in polarized hadron hadron collisions, High energy physics and quantum field 
theory, Tver (2000), p. 399, hep-ph/0107059| 

S. Dittmaier and M. Kramer, Phys. Rev. D65, 073007 (2002), |hep^ph /0109062[ 

F. Berends and R. Kleiss, Z. Phys. C27, 365 (1985). 

U. Baur and T. Stelzer, Phys. Rev. D61, 073007 (2000), |hep^h/9910206| 

E. Barberio, B. van Eijk, and Z. W§s, Comput. Phys. Commun. 66, 115 (1991), ibid. 
79, 291 (1994). 

S. Jadach, W. Placzek, M. Skrzypek, B. F. L. Ward, and Z. W§s, Comput. Phys. 



Commun. 140, 432 (2001), |hep-ph/0103163 



D. R. Yennie, S. Frautschi, and H. Suura, Ann. Phys. (NY) 13, 379 (1961). 

S. Jadach and W. Placzek, WINHAC 1.00: The Monte Carlo Event Generator for 
Single W-Boson Production with Leptonic Decays in Hadron Collisions, available 
from the authors; documentation in preparation. 

K. Hagiwara and D. Zeppenfeld, Nucl. Phys. B274, 1 (1986). 

D. Y. Bardin, A. Leike, T. Riemann, and M. Sachwitz, Phys. Lett. B206, 539 (1988). 

S. Haywood et al, Electroweak Physics, in 

J. Fleischer, F. Jegerlehner, and M. Zralek, Z. Phys. C42, 409 (1989). 

T. Kinoshita, J. Math. Phys. 3, 650 (1962). 

T. D. Lee and M. Nauenberg, Phys. Rev. B133, 1549 (1964). 

U. Baur and D. Zeppenfeld, Phys. Rev. Lett. 75, 1002 (1995). 

E. Argyres et al, Phys. Lett. B358, 339 (1995). 

M. Passera and A. Sirlin, Phys. Rev. D58, 113010 (1998), |hep^h/9804309| 

A. Sirlin, Recent developments in the study of unstable particle processes. Proceed- 
ings of the IVth International Symposium on Radiative Corrections (RADCOR 98): 
"Applications of Quantum Field Theory to Phenomenology", Barcelona, 1998, ed. 
by J. Sola (World Scientific, Singapore, 1999), p. 546, ,hep-ph/9811467, 



25 



[34] S. Jadach and B. F. L. Ward, Comput. Phys. Commun. 56, 351 (1990). 
[35] W. Placzek and S. Jadach, in preparation. 

[36] S. Jadach, W. Placzek, M. Skrzypek, and B. F. L. Ward, Phys. Rev. D54, 5434 
(1996). 



26 



